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mic parameters 

G.A. Verner* and I.W. Roxburgh 

Astronomy Unit, Queen Mary, University of London, Mile End Road, London, El 4NS, UK. 
The dates of receipt and acceptance should be inserted later 



o 

<N 
< 



go ; 
o : 

c« ■ 



m 

O 

o 




1000 2000 3000 4000 5000 6000 
Frequency (^Hz) 



Key words stars: oscillations - methods: data analysis 

This article describes a moving-windowed autocorrelation technique which, when applied to an asteroseismic Fourier 
power spectrum, can be used to automatically detect the frequency of maximum p-mode power, large and small sep- 
arations, mean p-mode linewidth, and constrain the stellar inclination angle and rotational splitting. The technique is 
illustrated using data from the CoRoT and Kepler space telescopes and tested using artificial data. 
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1 Introduction 

Regions of p-mode power can be characterised by the self- 
repeating comb-like structure in a Fourier power spectrum. 
Using a moving-windowed autocorrelation technique ap- 
plied to the acoustic power spectrum it is possible to au- 
tomatically detect the frequency of maximum power, large 
and small separations, mean mode linewidth and, with suf- 
ficient signal-to-noise, constrain the inclination angle and 
rotational splitting. With the rapidly increasing amount of 
asteroseismic data becoming available from the CoRoT and 
Kepler space telescopes, there is a need for accurate and 
flexible analysis pipelines which can automatically obtain 
global asteroseismic parameters. A number of automated 
pipelines have recently been described by groups preparing 
to analyse data from the Kepler mission (Hekker et al. 2010; 
Huber et al. 2009; Mathur et al. 2010). The autocorrelation 
technique described here has been applied to high and low 
signal-to-noise data from the CoRoT spacecraft and tested 
successfully with artificial data to show that it provides a 
powerful method to detect global p-mode parameters and is 
particularly useful for data with low signal-to-noise. 



2 Pipeline 

The first stage of the analysis pipeline is to account for 
the background variation in each power spectrum. A simple 
background estimate was found to be sufficient to remove 
the smooth variation in power without attenuating the sig- 
nal from the p-mode region. Therefore, the background was 
obtained from the power spectrum by applying a moving- 
median filter with a width of 100/iHz to obtain an estimate 
of the mean power, followed by a moving-mean filter with 
the same width to smooth the background estimate. The 
power spectrum is then scaled by this background to give 
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Fig. 1 Background fit and removal for HD49933. Top 
panel shows the background estimate (green) which has 
been removed in the scaled power-to-background spectrum 
below. 



a relative power-to-background spectrum. This is illustrated 
in Fig.[[]for the case of the CoRoT F2V star HD49933 (Ap- 
pourchaux et al. 2008; Benomar et al. 2009). 

Once the background has been removed, the windowed 
power- spectrum autocorrelation function (PS AC) is gener- 
ated from the relative power-to-background spectrum (PBS) . 
In each window of width 500/iHz (incremented in steps of 
1/iHz), the unbiased PS AC is calculated from (PBS - 1). 
Regions with no autocorrelation signal (i.e. no p-mode de- 
tection) return a PS AC of approximately zero. Those win- 
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Fig. 2 PSAC maps and z/ m ax fits for HD49933 (left), HD181906 (middle) and KIC6603624 (right). Top panels show 
the windowed PSAC as a function of central frequency, centre panels show the absolute values of the PSAC map, bottom 
panels show the mean absolute PSAC at each central frequency (black) and a Gaussian fit to the peak (red). 



dows containing p-mode signal show a positive and nega- 
tive correlation as the ridges alternate into and out of phase. 
By summing the absolute deviations of this PSAC map, the 
collapsed signal can be fitted to a Gaussian function to give 
the frequency of maximum power (^max) and the width of 
the p-mode region. This is illustrated in Fig. [2] in the case 
of HD49933, HD181906 - a low signal-to-noise F8 star ob- 
served by CoRoT (Garcia et al. 2009) and KIC6603624 - a 
G-type star in the first public release from Kepler (Chaplin 
etal. 2010). 

When the central p-mode region (z^max) has been iden- 
tified, an initial estimate of the large separation can be ob- 
tained from the narrow- windowed autocorrelation of the time 
series (TSAC) (Roxburgh 2009a). Further information can 
then be obtained from a fit to the PSAC in in a single win- 
dow containing the central region of maximum power. The 
amount of information available from the PSAC depends 
on the properties of the modes, and particularly the mode 
linewidths in relation to the small separations. The PSAC 
is shown in Fig. [3] for three artificial scenarios in different 
linewidth regimes. 

The PSAC can be fitted to the limit- spectrum PSAC ob- 
tained from a sum of Lorentzian peaks parametrised by the 
large separation (Az/), small separations (doi, ^02, ^13), av- 
erage mode linewidth, rotational splitting, inclination angle 
and height ratios. Some of these parameters can be fixed, or 
modes of degree t = 3 can be omitted, based on the signal- 
to-noise ratio and a priori information. The PSAC fits to five 
spectra are shown in Fig. [4] The fit to solar data (from the 
VIRGO/SPM instrument on the SoHO spacecraft) identifies 
the i = 3 peak and illustrates the narrow solar line profiles. 
The C0R0T F stars have large mode linewidth which limits 
the extraction of small frequency spacing. The Kepler G- 
type stars have narrow, solar-like linewidths and clear small 
frequency separations. In the case of KIC3656476, the cen- 
tral peak has a triple-peaked appearance due to the close 
similarity between the values of doi and <io2- 
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Fig. 3 Artificial normalised PSAC functions showing the 
doi separation (green dashed line), <io2 separation (blue 
dashed line) and Az/ separation (red dashed line). The ar- 
tificial spectra were generated with linewidths of 3.0, 1.0 
and 0.5/iHz (from top to bottom). The rotational splitting 
is non-negligible only in the bottom panel, where it is set 
at 0.5/iHz. The large separation was set at 100/iHz in each 
case. 

3 Testing 

The pipeline was tested with 177 artificial Kepler-like time 
series each of 30-day s duration (Mathur et al. 2009). The 
distributions of fitted parameters about their input values 
are shown in Fig. [5] The determination of vmax returned a 
value in 83% of cases. The cases where no value was found 
corresponded to artificial data with amplitudes of less than 
3ppm. The accuracy of the returned value of z/ ma x shows 
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Fig. 4 Fits to the PS AC function centred on ^max for 
five spectra - Sun (VIRGO/SPM), HD49933, HD181906, 
KIC6603624 and KIC3656476 (from top to bottom). Red 
lines show the fit to the data. 



a dependence on the amplitude of the modes. For ampli- 
tudes up to 4ppm, the mean absolute deviation from the true 
value was ~ 100/iHz, but this reduces to ~ 40/iHz for am- 
plitudes of 8ppm and less than 10/iHz at lOppm. The large 
spacing determination from fitting the PSAC shows an im- 
provement in the accuracy of the value from that obtained 
from the TSAC. Again, there is a dependence on amplitude, 
but for the lowest amplitude data, the mean absolute de- 
viation in fitted Au decreases from ~ 4/iHz in the TSAC 
case to ~ 2/iHz in the PSAC case. For the small frequency 
spacings (<ioi and ^02) the accuracy is more strongly de- 
pendent on the mean 'observed' mode linewidth than the 
maximum mode amplitude. The 'observed' mode linewidth 
is a combination of the linewidths of individual multipole 
components (due to damping and excitation) and rotational 
splitting. The doi spacing was found to be more readily ob- 
tained and more accurate than <io2- In cases where the mean 
linewidth of the modes was less than 1/iHz, the mean ab- 
solute deviation in fitted doi values was less than 0.4/iHz, 
compared with 1.2/iHz for <io2- 
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Fig. 5 Results of fits to artificial data. Left panels show the 
scatter of the fitted values from their real values, right panels 
show the histogram of the absolute deviations from the true 
value for each parameter with an error bar indicating the 
standard deviation in each bin. Results are shown for the 
fits to (from top to bottom): ^max, Av (from TSAC fit), Az/ 
(from PSAC fit), d§\ and <io2 small separations. 



www.an-journal.org 



© WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



4 



G.A. Verner & LW. Roxburgh: Power spectrum autocorrelation techniques to detect asteroseismic parameters 



4 Conclusions 

The power- spectrum autocorrelation provides a useful way 
to reliably detect global asteroseismic parameters in an au- 
tomated way. For the best candidates, a full power- spectrum 
maximum-likelihood fit obtaining individual mode parame- 
ters is preferable, but in cases where this is not possible, or 
is complicated by low signal-to-noise, obtaining an accurate 
estimate of the large and small separations is desirable. For 
cases where full power- spectrum fitting is possible, it is still 
necessary to generate initial guess parameters for the fitting 
model, which can be estimated more easily if the average 
global parameters are known. 

Using the PS AC technique, it is found that the doi small 
separation is more easily determined than the more common 
do2 small separation due to the influence of the doi spacing 
creating a distinctive double or triple peak in the PS AC ob- 
tained at ^max- It has been shown that the doi spacing is 
sensitive to the structure of the inner core and convective 
envelope (Roxburgh 2009b), it is therefore reassuring that 
this parameter may be detectable in cases where it is not 
possible to obtain a robust estimate of the <io2 spacing. 
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